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(PRELIMINARY) 

MODELING, SIMULATION AND CONTROL 
FOR 

A CRYOGENIC FLUID MANAGEMENT FACILITY 


How relevant to the complex space station applications are academic 
formulation and solution of control problems, based on recently published 
textbook methodology complemented with limited laboratory scale experiments? 
The textbook abstractions are often stripped of consideration of constraints 
of prime concern to the field application: process capacities, user demands, 

economics, hazards analysis and fault tolerance. However, the approach of 
the classroom — simplistic of necessity due to man-hour and funding 
constraints — serves as a starting point for formulating a "top-down 
modular" definition of the problem and development of an overall perspective 
for the research professor or student. The individual is thus conditioned to 
readily adapt to a position in team efforts with major funding. 

As one of an ongoing series of term projects in Process Monitoring and 
Control at UH-CL, the class in PROC 5232: Process Modeling, Simulation and 
Control, has studied the synthesis of a control system for a cryogenic fluid 
management facility. The severe demands for reliability as well as 
instrumentation and control unique to the space station environment are prime 
considerations . 

Realizing that the effective control system depends heavily on 
quantitative description of the facility dynamics, a methodology for process 
identification and parameter estimation is postulated. A block diagram of 
the associated control system is also postulated. Finally, an on-line 
adaptive control strategy is developed utilizing optimization of the velocity 
form control parameters — - proportional gains, Integration and derivative 
time constants — in appropriate difference equations for direct digital 
control. 

Of special concern are the communications, software and hardware 
supporting interaction between the ground and orbital systems. It is 
visualized that specialists in the OSI/ISO utilizing the Ada programming 
language will influence further development, testing and validation of the 
simplistic models here presented for adaptation to the actual flight 
environment. 
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1 . PROCESS DESCRIPTION 


1.1 Baseline Configuration: 

CFMFE Flight System 

The initial concept is diagrammed in Figure 1.1. Assembled 
as a module for mounting in the shuttle, it consists of three 
submodules identified with successive operational stages: 

a) chilldown of the Ground Fill Line on the pad; 

b) chilldown and filling of the Supply Tank on the pad; 

c) chilldown of the Transfer Line combined with chilldown 
and filling of the Receiver Tank in orbit. 


The submodules for operational stages a) and b) are detailed in 
Figure 1.1.1. The submodules for operational stages c) and d) 
are detailed in Figure 1.1.2. 


NOTICE 


At the deadline for submitting manuscripts this paper was 
incomplete . 

Copies of the completed version will be made available at 
the presentation to those who desire one. 


Max Turner 

UH-Clear Lake Box 329 
2700 Bay Area Blvd. 
Houston, TX 77058-1098 
(713) 488-9480 
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FIGURE 1.1.1 GROUND FILL OF LH2 SUPPLY TANK 
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1. PROCESS DESCRIPTION 

1.2 Functional Requirements and Constraints 


Time: 

Chilldown time. 

e = 

c 

15 min. 


Fill time. 

0 F ” 

60 min. 

Pressure: 

Minimum 

P min 

■ 1 atm 


Maximum excursion 

P 

max 

■ 85 psig (PSV spec) 

Temperature: 

Minimum 

T 4 

min 

*= 36.7°R (20. 4°K) 


Ambient 

T = 
a 

530°R 


Maximum ATg^ 

AT 87max " < TBD > 


Conservation of H2: (TBD) 


Hazards: 


T 


Explosion and fire (TBD) 
Destructive vibration (TBD) 
and shock 

Stress fractures (TBD) 

Loss of power (TBD) 


Zero-gravity: Liquid pressurization (TBD) 

Chilldown of receiver tank system (TBD) 

Filling receiver tank (TBD) 

Contingency respondent and fault tolerant (TBD) 
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1. PROCESS DESCRIPTION 

1.3 Problem Identification 

1.3.1 Thermal balances and minimized system chilldown and fill times 
On the pad: 

1. The Ground Fill line 

2. The CH2 Storage and Supply Tank 

In orbit: 

3. Chilldown time for the transfer line from the Supply Tank to the 
Receiver Tank 

4. Chilldown and fill time for the Receiver Tank 
Special problems: 

5. Overpressures and destructive stresses 

6. Delayed GH2 boiloff due to heat transfer 
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2. MATHEMATICAL MODELING 

2.2 Analysis, Degrees of Freedom and Control Loops 

2.2.1 The Ground Fill Line 

Physical model for an energy balance: 

Assumptions : 

1. The aluminum tube is perfectly insulated 

2. LH2 enters with quality x - 0 

3. Until chilldown is essentially complete, the exit GH2 has a 
quality of x - 1 

4. Maximum chilldown rates are limited by the venting capacity of the 
line 

5. Significant thermal energy sources which limit the minimum 
cooldown time are the concentrated masses associated with 
stainless steel control valves and sensors. 

6. The enthalpy of LH2 at near atmospheric pressure is given by: 

h = 278.4 + 441. 8x + 10.13 (T-21) 

■ 507.47 + 10.13T, kJ/kg using the unit °K 

- 218.63 + 2.425T, Btu/lb using the unit °R 

in 

Reference: Perry and Green, Ch.E. Handbook, McGraw-Hill 1984, 

pp 3-1958 

C , - C -10.13 kJ/kg°K at 21°K 
Pf Pg e 

7. The heat capacity of A1 is: 

C vA1 - - 0.1362 + 0.007528T - 0.00001356T 2 kJ/kg°K 

with T in °K 

C A . - - 0.03254 + 0.000999T - 9.99 x 10" 7 T 2 Btu/lb °R 
vAl m 

with T in °R 

Reference: Perry and Green, 1984, pp 3-135 

8. The heat capacity of stainless steel is: 

C v - - 0.0586 + 0.003219T - 5.078 x 10~ 6 T 2 kJ/kg°K 

with T in °K 

C - - 0.0140 + 0.000428T - 3.75 x 10“ ? T 2 Btu/lb °R 
v m 

with T in °R using 1 Btu/lb^R - 4.178 kJ/kg°K 
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2.2.1 The Ground Fill Line 


Assumptions (continued) 


9. Radiation heat transfer rates across the annulus of concentric 

tubes or spheres is nil compared to convective heat transfer rates 
from A1 to LH2: 

q/A = 300 Btu/hrft 8 from A1 to LH2 at 36.7°R 


q/A 


12 


e r (T - T LH2 ) 


12.7 Btu/ft hr from StSt to A1 


at 530 R where F, 


/ 


/ 


12 

& - .1713 x 10 8 Btu/(f t^hr°R^) 


7 


£ 6 . 69 ^ 


Reference: Perry and Green, 1984 


10 . 


Thermal dif fusivities , k/pC v : 

Al 

St. St. 

20°K 

0.5 

0.040 

100°K 

0.00023 

~ 9 x 10 -5 

300°K 

0.00011 

3.3 x 10 


Reference: Perry and Green, 1984, pp 3-263 

/ /, ^GH2 ° °* 10 ^ kg/® 3 at 20.4°K 

p TH ? = 70.57 kg/m 8 at 20.4°K or 4.72 16 m /ft 8 

3 3 

11. Densities: p A1 *= 2723 kg/m or 170 lb /ft 

' Ai nx 

pg t gt “ 7900 kg/m 8 or 492 lb m /ft 8 
Reference: Perry and Green, 1984, pp 3-96 

12. Thermal conductivities: 

Reference: Perry and Green, 1984, pp 3-261 

13. Convective heat transfer coefficients: 


References: 1) Perry and Green, 1984, pp 10-23 

2) H.H. Walters, AiResearch Manufacturing Compant 
"Single-Tube Heat Transfer Tests with Liquid 

Hydrogen", (see WADC Technical Report 59-423) 

3) Drake et al. , Arthur D. Little, Inc. 
"Pressurized Cool-Down of a Cryogenic Liquid 

Transfer system Containing Vertical Sections", 
(test£ yi|h L02) 
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2.2.1 The Ground Fill Line 


Assumptions (continued) 

Walters — LH2 tests: 

film boiling: h D 460 to 540 Btu/hrft^°R 

for inlet (?) = 1.6 to 1.7 atm 

Re • 3 x 10 5 

nucleate boiling: h = 10 x value for film boiling 

Drake et al. — L02 tests: 

film boiling: h = 300 Btu/ft^hr°R 

for inlet pressure = 20 psig 
outlet pressure = 10 psig 

h = 200 Btu/ft^hr°F 
for inlet pressure “ 10 psig 
outlet pressure =5.5 psig 

LH2 — - assume: h = 500 Btu/hrft^°R (uncontrolled) 

max 

h = 300 Btu/hrft^°R (controlled) 
ave 

Re = 3 x 10 5 

14. Critical constants of H2: P » 12.8 atm 

c 

T = 33.3°K = 60.0°R 
c 
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